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Abstract 

Trapped dynamics widely appears in nature, e.g., the motion of particles in viscous cyto¬ 
plasm. The famous continuous time random walk (CTRW) model with power law waiting 
time distribution ( having diverging first moment ) describes this phenomenon. Because of 
the finite lifetime of biological particles, sometimes it is necessary to temper the power law 
measure such that the waiting time measure has convergent first moment. Then the time 
operator of the Fokker-Planck equation corresponding to the CTRW model with tempered 
waiting time measure is the so-called tempered fractional derivative. This paper focus on 
discussing the properties of the time tempered fractional derivative, and studying the well- 
posedness and the Jacobi-predictor-corrector algorithm for the tempered fractional ordinary 
differential equation. By adjusting the parameter of the proposed algorithm, any desired 
convergence order can be obtained and the computational cost linearly increases with time. 
And the effectiveness of the algorithm is numerically confirmed. 
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1 Introduction 


The fractional calculus has a long history. The origin of fractional calculus can be 
traced back to the letter between Leibniz and L’Hopital in 1695. In the past three 
centuries, the development of the theories of fractional calculus is well contributed 
by many mathematicians and physicists. And from the last century, the books cov¬ 
ering fractional calculus began to emerge, such as Oldham and Spanier (1974), 
Samko, Kilbas and Marichev (1993), Podlubny (1999), and so on. In recent years, 
more theories and experiments show that a broad range of non-classical phenomena 
appeared in the applied sciences and engineering can be described by fractional cal¬ 
culus H33I26I35L Because of its good mathematical features, nowadays fractional 
calculus has become a powerful tool in depicting the anomalous kinetics which 
arises in physics, chemistry, biology, finance, and other complex dynamics If26tl . 
In practical applications, several different kinds of fractional derivatives, such as 
Riemann-Liouville fractional derivative, Caputo fractional derivative 11331351 . Riesz 
fractional derivative ll35l . and Hilfer fractional derivative H191421 are introduced. 

One of the typical applications for fractional calculus is the description of anoma¬ 
lous diffusion behavior of living particles; and the tempered fractional calculus de¬ 
scribes the transition between normal and anomalous diffusions (or the anomalous 
diffusion in finite time or bounded physical space). In the continuous time random 
walk (CTRW) model, for a Levy flight particle, the scaling limit of the CTRW with 
a jump distribution function (p(x) ~ x~ (A+<,) ( \ < a < 2) exhibits superdiffusive dy¬ 
namics. The corresponding stable Levy distribution for particle displacement con¬ 
tains arbitrarily large jumps and has divergent spatial moments. However, the infi¬ 
nite spatial moments may not be feasible for some physical processes [[8J. One way 
to overcome the divergence of the moments of Levy distributions in transport mod¬ 
els is to exponentially temper the Levy measure. Then the space fractional operator 
will be replaced by the spatially tempered fractional operator in the corresponding 
models mm. This paper concentrates on the time tempered fractional deriva¬ 
tive, which arises in the Fokker-Planck equation corresponding to the CTRW model 
with tempered power law waiting time distribution B4I17H . Tempering the power 
law waiting time measure makes its first moment finite and the trapped dynamics 
more physical. Sometimes it is necessary/reasonable to make the first moment of 
the waiting time measure finite, e.g., the biological particles moving in viscous cy¬ 
toplasm and displaying trapped dynamical behavior just have finite lifetime. The 
time tempered diffusion dynamics describes the coexistence/transition of subdif¬ 
fusion and normal diffusion phenomenon (or the subdiffusion in finite time) which 
was empirically confirmed in a number of systems 1181271 . More applications for the 
tempered fractional derivatives and tempered differential equations can be found, 
for instance, in poroelasticity [H8ll . finance Q, ground water hydrology Il27l28l . 
and geophysical flows ||291 . 

Tempered fractional calculus can be recognized as the generalization of fractional 
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calculus. To the best of our knowledge, the definitions of fractional integration with 
weak singular and exponential kernels were firstly reported in Buschman’s earlier 
work a. For the other different definitions of the tempered fractional integration, 
see the books H39I35I281 and references therein. This work continues previous ef¬ 
forts lf25l to explore the time tempered fractional derivative. The well-posedness, 
including existence, uniqueness, and stability, of the tempered fractional ordinary 
differential equation (ODE) is discussed, and the properties of the time tempered 
fractional derivative are analyzed. Then the Jacobi-predictor-corrector algorithm 
for the tempered fractional ODE is provided, which has the striking benefits: 1. any 
desired convergence order can be obtained by simply adjusting the parameter (the 
number of interpolation points); the computational cost increases linearly with the 
time t instead of t 2 usually taken place for nonlocal time dependent problem. And 
extensive numerical experiments are performed to confirm these advantages. 

In Section [2} we introduce the definitions and show the properties of the tempered 
fractional calculus, including the generalizations of the tempered fractional deriva¬ 
tives in the Riemann-Liouville and Caputo sense, and the composite property. More 
basic properties are listed and proved in Appendix A; the expression and proper¬ 
ties of the tempered fractional calculus in Laplace space are proposed and proved 
in Appendix B. In Section [3j we discuss the initial value problem of the tempered 
fractional ODE: first derive the Volterra integral formulation of the tempered frac¬ 
tional ODE; then prove the well-posedness of the considered problem. The Jacobi- 
Predictor-Corrector algorithm for the tempered fractional ODE is designed and dis¬ 
cussed in Section |4] and two numerical examples are solved by the algorithm to 
show its powerfulness. 


2 Preliminaries 

In this section, we first give the definitions and some properties of the tempered 
fractional calculus. Let [a, b ] be a finite interval on the real line R. Denote L([a, b\) 
as the integrable space which includes the Lebesgue measurable functions on the 
finite interval [a, b ], i.e., 



And let AC[a, b ] be the space of real-values functions u(t) which are absolutely 
continuous on [a,b\. For n e N + , we denote AC"[a, b] as the space of real-values 
functions u{t) which have continuous derivatives up to order n — 1 on [a, b\ such 
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And denote by C"[a, b ] the space of functions u{t) which are n times continuously 
differentiable on [a, b\. 


Definition 1 (Riemann-Liouville tempered fractional integral ||4f8l ) Suppose that the 
real function u{t) is piecewise continuous on ( a , b ) and u{t) £ L([a, b ]), cr > 0, A > 

0. The Riemann-Liouville tempered fractional integral of order cr is defined to be 

J? A u(t) = e~ At a I?(e A, u(t)) = —— f e~ A( ‘~ s) (t - sf~ ] u(s)ds, ( 1 ) 

r(cr) J a 

where a lf denotes the Riemann-Liouville fractional integral 


J>(t) 


=f hlfi 


■S')' 7 ' l u(s)ds. 


( 2 ) 


Obviously, the tempered fractional integral ([[]) reduces to the Riemann-Liouville 
fractional integral if A = 0. In practical applications, sometimes the fractional inte¬ 
gral (□} is represented as a D~ rT ’ A u(t). 

Definition 2 (Riemann-Liouville tempered fractional derivative BI8I0 For n - 1 < 
a < n,n e N + , A > 0. The Riemann-Liouville tempered fractional derivative is 
defined by 


a D a Mt ) = e- M a D a t (e A ‘u(t)) 


At d n r t 


f 


e u(s) 


-ds, 


T{n - a) dt n J a (t - s) a " +1 
where a D'' (e At u(t)) denotes the Riemann-Liouville fractional derivative f\33\] 


( 3 ) 


d n 


M;(e lt u(t)) = —c/rvw/))) = 


1 d n 


n nt t„As 


f 


( e As u(s )) 


Y(n - a) dt n J a {t - s) a ~ n+ 


-ds. (4) 


Remark 1 fll3ln The variants of the Riemann-Liouville tempered fractioncd deriva¬ 
tives are defined as 


Df A u(t) = 


a Df A u(t) - A a u(t), 

,,Df A u(t) - - A a u(t), 

dt 


0 < a < 1, 
1 < a < 2 . 


( 5 ) 


Definition 3 (fractional substantial derivative Ill6l40l6l0 For n —1 < a < n,n e N + , 
and A(x) being any given function defined in space domain. The Riemann-Liouville 
fractional substantial derivative is defined by 


= (4 + «) a ir M u(t) = (4 + aw) f (6) 
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where ,r " Mx] denotes the Riemann-Liouville fractional integral and 


j, +m ) 


(V) 


Remark 2 The fractional substantial derivative d6]) is equivalent to the Riemann- 
Liouville tempered fractional derivative © if A(x) is a nonnegative constant func¬ 
tion. In fact, using integration by parts leads to 


(I + m 
= (I + 4j:) 

= {jp^) 

= (| + AM 


f 


■t g -A W-(/-i) M ( s ) 




« a - 5)“-" +1 

n ~ lrt d 

- 1 - A(x). . 

dt ) J a (t-s ) a ~ n+1 


t e -A(x)-(t-s) u ( s ) 


ds 


n— 1 


n -2 


d C e Mx)s u(s) 


e~ A(x)t - I -- '^1—ds 

(t - s) a ~ n+l 


-A(x)t 


d_ r 

dt J a 

<fi_ f -J— 

dt 2 J a (t - 


A e A(x)s u{s) 


sy 


,a-n+\ 


ds 


( , \n 

j t + A\ u(t), which can be 

simply/resonably denoted as D"- A u(t). 


Definition 4 (Caputo tempered fractional derivative Ii35l41l0 Forn- 1 < a < n,n e 
N + , A > 0. The Caputo tempered fractional derivative is defined as 


Cr> a ’ A 

a 


D^u(t) = e~ At c a D a t {e At u(t )) = 


„-At 


Y(n - a) 


r 1 

J a (t-sr~" +i 


d n (e As u(s)) 


ds n 


ds, (8) 


where c a D‘J A (e A, u(t )) denotes the Caputo fractional derivative 1 13311 


]D a Ae At u{t )) =--- f ---- 

r (n-a)J a (I - s f~ n+l 


d n (e As u(s)) 


ds n 


ds. 


(9) 


Remark 3 The equivalent forms of Riemann-Liouville tempered fractional deriva¬ 
tive (0) and Caputo tempered fractional derivative © are a D°f ,A u(t ) = D"- A a If a ' A u(t) 
and c a D a f X u{t) = a If~ a ’ A D n,A u{t), respectively. 

Note that when A = 0, the Riemann-Liouville (Caputo) tempered fractional deriva¬ 
tive reduces to the Riemann-Liouville (Caputo) fractional derivative. 


Proposition 1 Let u{t) e AC n [a, b\ and n - 1 < a < n. Then for all t > a, holds 


c a Df\u{t )) = a Df\u{t)) - Y, 


— e~ A, (t - a) k ~° d k 


k=0 


T(k-a+ 1) 




( 10 ) 
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Proof. Take v(t) = e Al u(t) in the equation for the Riemann-Liouville and Caputo 
fractional derivatives It33i23l35fl 


c 

a 


D a t v{t) = a D a t Ut) - Yj 


n— 1 


k=0 


(t - a) k d k v(t) | \ 

~~k\ 


yielding 


77—1 


D L f(e At u(t)) = a D°[e At u(t) - Y ^ ^ “(0)L,)- 


k j k 


k =0 


Multiplying both sides of the above equation by e , we obtain 


77—1 


e A ‘ a D “(e A, u(t )) = e At a D°[e At u(t) - Y ^ 


k=0 


Furthermore, using the definitions of Riemann-Liouville and Caputo tempered frac¬ 
tional derivatives, we get that 


c 

a 


Df\u(t)) = a Df\u(t)) - 


a. A/ 



(ll) 


Using the linearity properties presented in Proposition [4] and the formula of power 
function 



(t ~ a) k \ 

k\ I 


i a D a t ({t-a) k 


1 T(k + 1)(7 - a) k ~ a 
k\ r(k - a + 1) 


we deduce the desired result from (fill) . 


□ 


Proposition 2 (Composite properties) 

(1) Let u(x) 6 L([a,b]) and I n ~ a ’ A u(t ) 6 AC n [a,b]. Then the Riemann-Liouville 
tempered fractional derivative and integral have the composite properties 

77-1 

a If A [ a Df A u(t)] = u(t) ~Y~ 

k =o 


ft - a) 


a-k -1 


r(a - k) 


-[ a D“- k -\e At U (t))\ t=a l (12) 


and 

a Df A [ a If A u(t )] = u(t). (13) 

(2) Let u(t) € AC" [a, b\ and n — 1 < a < n. Then the Caputo tempered fractional 
derivative and the Riemann-Liouville tempered fractional integral have the com- 
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posite properties 


and 


n— 1 




at 


u{t) - Y 


k=0 


e 


A,H-a) k 

\d k (e A, u(t)) 

- 

k\ 

dt k 

t-a- 


c a D a t \l°Mm = u{t) if a e (0,1). 


(14) 


(15) 


Proof. From the definitions of Riemann-Liouville tempered fractional integral and 
derivative, we have 


r? A [ a Dfu(t)] = e- At a r;[e A \ a D a t Mt))] 


a L t 


= e - M a I?[e M (e- M a D?( e At u(t)))\ 
= e~ At X[ a D L ;{e M u{t))\. 


(16) 


(i) 


Thanks to the composition formula H33I231351 


n -1 


a IXD?v(t)] = v{t) - Y 


(t - a) 


a—k—l 


k =0 


r(a - k) 




we get 


(7) = a r*[ a D«e M u(t)} = e M u(t) - Y 


(t - a) 


a-k-l 


k =0 


r(a - k) 


-[ a Dr k ~\e At 


“W)| tJ- 


Inserting the above formula into (fl6l) leads to (fl2l) . 


Again from the definitions of Riemann-Liouville tempered fractional integral and 
derivative, there exists 


a 




= e- At a D?[e A ‘U?’ A um 

= e~ At a D" [e At (e~ At a If (e At u(t)))\ 

= e- At a D c ;[X(e At u(t))]- 


Furthermore, using the composite properties of Riemann-Liouville fractional inte¬ 
gral and derivative B3l23i35l 


a D a t [JtW))] = HO, 


(17) 


we get 

aDflXMm = e~ A, Hy;\JHe A, u(0)\ = U{t), 
by taking v(t) = e A, u{t) in Eq. (fT71) . 
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Similarly, using the composite properties of Caputo fractional derivative H33I23I35II 


n— 1 


= v(f) - £ 

k= 0 


(r-a)* r rf*v(Q i 

fc! 1- dt k t=a 


and 

^Z)“[ a /“(v(0)] = v(f), if a € (0,1), (18) 

we can get (fT4l) and (fl5l) . 


□ 


Remark 4 For a constant C, 

a D a /C = Ce~ At a D c ;e At , c a Df A C = Ce~ M c a D“e M . (19) 

Obviously, a D n l ~ A (C.) =£ ^D"' A (C). And is no longer equal to zero, being 

different from ( u D"(C) = 0. 


3 Well-posedness of the tempered fractional ordinary differential equations 


In this section, we consider the ODEs with Riemann-Liouville and Caputo tem¬ 
pered fractional derivatives, respectively, i.e., 


a Df A u(t) = f{t, u(t)), n - 1 < a < n, A > 0, 

[ a Dff k ~\e At um\ t=a = 8k, k = 0 , 1 , 2 , • • • ,n - 1 , 


( 20 ) 


and 


c a D° t ' u(t) = fit, u(t)), n - l < a < n, d>0, 
d k 1 

—( e At u(t )) = c k , k = 0,1,2, ■ • • ,n - 1. 

Clt t-a 


( 21 ) 


The Cauchy problems (l20l) and (12TI) can be converted to the equivalent Volterra 
integral equations of the second kind under some proper conditions. 


Lemma 1 If the function f(t, u(t)) and u{t) belong to L([a, b ]), then u{t ) is solution 
of the initial value problem (l20l) if and only if u(t) is the solution of the Volterra 
integral equation of the second kind 


u(t) = ^ 


— e~*(t - a) a ~ k - 1 


8 k- 


+ 


k=0 


T(a - k) r(or) 


l r 

Ja 


-A(t-s)/ f _ p'vQ'-l 


(t - s) a fis, uis))ds. (22) 


In particular, if 0 < a < 1, then uit ) satisfies the Cauchy problem (1201) if and only 
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if uit) satisfies the following integral equation 


n —At 


U{t) = ga¬ 


it-a) 


T(a) 


a -1 2 

+ Ha) 1, ' 


-A(t-s) 


it 


sf l fis,uis))ds. 


(23) 


Proof. For the linear Cauchy problems ([20b and (l2TT) . the conclusion is directly 
reached by the Laplace transform given in Appendix B. Now we prove the more 
general case. 

Necessity. Performing the integral operator a If A on both sides of the first equation 
of (l20l) . we have 


n —1 


«(0 = Y 


-At 


gk- 


it - a) 


k=0 


r ia - k) 


a-k-\ ^ nt 


-Mf-s) 


it 


sf l fis,uis))ds. 


where we use the composite property (1) given in Proposition [2] Then Eq. (l22l) is 
obtained. 


S ufficiency. Applying the operator a D‘f A to both sides of Eq. (l22l) results in 


n -1 


uDfMt) = Yj Sk — {( \, , (/ , Uf k l) + aDf A a If A fit, Uit )) = fit , uit)), (24) 




k =0 


T(a - k) 


where we use the fact 

a Df A ie~ At it-af- k - 1 ) e~ Al it - af- k - l ~ a e~ A, it - a)~ k ~ l 


rtcr - k) 


r i-k) 


oo 


0, k = 0,1,2, • • • , 77-1, 


and the composite property (fl3l) . Now we show that the solution of (l22l) satisfies the 
initial condition given in Eq. (l20l) . Multiplying e At and then performing the operator 
on both sides of Eq. (1221) . for 0 < j < n — 2</7 — 1 < cr < n, we have 


n —2 


B Z>r - V k(0) = X - . — + aD c r H afie A, fit, uit))). 


k =0 
n —2 


=z 

k=0 


gk 


Tij-k+l) 

it - af~ k 
Tij-k+l) 


(25) 


+ fl D^- 1 (e^/(r,«(0)), 


where the formula 


a D°f A ~ 1 (((- af- k ~ 1 ) 


//“L 1 


iit-af-*- 1 ) 


Via - k) 
Tij-k+l) 


it-af k , 0 < k < 77-2, 


is utilized. 
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Taking a limit t —» a in the above equation, we obtain 

«-2 , \j-k 

\im a D a t ~ ] ~\e At u{t)) = lim V g k ^ . ° — + lim a D~ ] ~\e At fit, u(t))), (26) 

t-*a t-*a 1 ( ] — k + 1) f-»a 

with the second term in the right hand side being equal to zero; and for its first term, 
we have 


n —2 


lim V 

t->a 


gk- 


(.t - ay 


-k 


./■-l 


f-a ^ T (j -k + 1) 

Jt=0 w ' 


= hm > g* ——r — + gj + lim > g* 
t—*o T(j - k + 1) g—i 


k =o 


^ r( j-k + i) 

k=j +1 U ' 


7-1 


gk 


n—2 


k =0 

= gj 


-inj-k+D 


o + gj + lim V 


gk(t - ay 


i-k 


k=j +1 


(27) 


□ 

By the similar technique in proving Lemma [Q we obtain the following conclusion 
for the Cauchy problem (OTI) . 


Lemma 2 If the function f(t, u) is continuous, then u{t) is the solution of the initial 
value problem (l2TT) if and only ifujt ) is the solution of the Volterra integral equation 
of the second kind 


u(t) = y 


— e At (t - a) k 


c k - 


+ 


k =0 


T(k + 1) T(a) 


1 f e- A(t -yt-s) a -'f(s,u(s))ds. (28) 
(«) Ja 


In particular, if 0 < a < 1, then u{t) satisfies the Cauchy problem if and only ifu{t) 
satisfies the following integral equation 


u(f) = u{a)e~ A(t - a) + — f e~ A(, - s \t - 
Ha?) Ja 


s) a f(s,u(s))ds. 


(29) 


3.1 Existence and uniqueness 


Many authors have considered the existence and uniqueness of the solutions to the 
nonlinear ODEs with fractional derivatives; see, e.g., H32I1I15I14I13I20I211421451 . 
For the global existence and uniform asymptotic stability results of fractional func¬ 
tional differential equations corresponding to (1231) . one can see [25121 • In the follow¬ 
ing, we discuss the existence and uniqueness of the solutions of the nonlinear tem¬ 
pered fractional differential equations based on the equivalent Volterra equations 
presented in Lemmas Q] and [2j We shall employ the Banach fixed point theorem to 
prove it. Let / : [a, b] xB —» R be a continuous function such that fit, u) e L([a, b ]) 
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for any u e B, being an open set in R. In the following, we always suppose that 
fit, u) satisfies the Lipschitz type condition with respect to the second variable 

I fit, u) - fit , v)| < C Lip \u - v|, for all u,v G B, t G [a, b ], (30) 

where C Lip is constant. We shall use the following space 

L a \[a,b ]) = ju G Li[a,b]), a Df A up) € IX[a,b])\. 


Theorem 1 Ifn - 1 < a < n,n G N + , A > 0, then there exists a unique solution uit) 
to the Cauchy problem (l20l) in the space L a ’ A i[a, b\). 


Proof. The proof of this theorem is similar to the references H32I14120142H . First we 
prove the existence of a unique solution uit) G L([a, b ]). In accordance with Lemma 
Ul it is sufficient to prove the existence of a unique solution uit) G L([a, b]) to the 
nonlinear Volterra integral equation (l22l) . We rewrite the integral equation (1221) in 
the form of operator 

uit) = (Pu)(t ), 


where 


and 


iPu)it) = uoit) + —— f e A(t s \t - sf l fis, uis))ds, 

n<*) Ja 


n —1 

Uoit) = 2 8k 


k =0 


e~ Al it - df ~ k - 1 
Via - k) 


(31) 


(32) 


We first prove that P is a contraction operator in the subinterval [a, t\\ c [<s, b] ia < 
t\ < b). To do this, we select t\ G ([o, b ]) such that the inequality 


C 


Lip 


ih - af 
Via + 1 ) 


< 1 


(33) 


holds.To apply the Banach fixed point theorem in the complete metric space Li[a, t\\), 
we have to prove the following facts: 


(i) If uit) G Li[a, ]), then (Pw)(t) G Li[a, t x ]); 


(ii) For all u\, u 2 g L(\a, /, J) the following inequality holds 


\\PUl - PUiWuaM) ^ ^Flll w l - U2\\L([aMb ^1 “ Clip + 1) ' (34) 

In fact, since fit, up)) G L([a, /,]) and Lemma |6] in Appendix A, the integral in 
the right-hand side of (PTT) belongs to L([a, q]); obviously u 0 p) e Li[a, q]), hence 
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(. Pu){t ) € L([a, /]]). Now, we prove the estimate (l34l) . From Lemma[6]in Appendix 
A, we obtain 

\\PUi - Pu 2 \\L«a M ) =11 air A fit, Mi(f)) ~ U 2 mWa M ) 

= U\(t)) ~ f(t, U 2 (t))) |li([a,fi]) 

<Wi||«i(0 - U 2 {t)\\ L ([a,h])- 

In view of 0 < W\ < 1 there exists a unique solution u*{t) £ L([a, ?!]) to Eq. (13TT) on 
the interval [a,t\\. The solution u(t) is obtained by taking the limit of convergent 
sequence ( P' n u* 0 )(t ) as m —» oo, i.e., 

lim || P m u* 0 - u*\\ mM) = 0, (35) 


where u* 0 (t) e L(\a, b\). Now let us consider the interval [q, t 2 \ with t 2 = t\ + h and 
h = t\ — a. Rewrite Eq. (1221) in the form 


u(t ) =u 0 (t ) + 


j_ r e -Mis ){t _ u(s))ds 

r(a) Ja 


+ ^~ f e- Mt - s \t-sT-\f(s,u(s))ds. 

r(Q') J tl 


(36) 


Since the function u(t) is uniquely defined on the interval \a,t{\, the last integral 
can be considered as the known function. Then the above equation can be rewritten 
as 


U(t) = ll()\ (t) + 


1 

r Xa) 



e W l f(s,u(s))ds. 


(37) 


«oi(0 = «o(0 + tt— e X(t s) (t - s) a l f(s, u(s))ds, (38) 

r(«) Jo 

is the known function. With the same contraction factor W\, we can prove that 
there exists a unique solution u*(t ) e L{t\,t 2 ) to Eq. (l22l) on the interval [fi, t 2 ]. By 
repeating this process finite times, e.g., M times, we can cover the whole interval 
[a,b]. 


To complete the proof of the theorem we must show that such a unique solu¬ 
tion u(t) £ L([a,b ]) belongs to the space L a ' A ([a, b]). It is sufficient to prove that 
a D a t ’ A {u(tj) £ L{[a, b ]). By the above proof, the solution u(t) £ L([a, b ]) is a limit of 
the sequence u m (t) £ L([a,b ]), i.e., 


lim || u m - u\\ U[aM) = 0, (39) 

m—>oo 

with the choice of certain u m on each [a, t{\, [fj, t 2 ], ■ • • , [t M -i,b]. In view of 
Ila77f U m aDf M||i([a,i>]) — \\f Ji U m -\) ~ f (t, U)\\L([a,bY) — CLip\\u m —l W||z.([a,£>]) , (40) 
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taking the limit of (1401) as m —> oo, gives 



(41) 


m 


oo 


and hence a Df A u(t) e L([a, b\). This completes the proof of the theorem. 


□ 


By almost the same idea, we can prove the following existence and uniqueness 
result for the Cauchy type problem (I2TI) . 

Theorem 2 If n - 1 < a < n, n s N + , A > 0, then there exists a unique solution 
u(t) to the Cauchy type problem (I2TI) in the space AC n [a, b\. 


3.2 Stability 

To prove the stability of the solutions of the Cauchy type problems (l20l) and (12T1) . 
we need the following generalized Gronwall’s Lemmas. 

Lemma 3 dl9H) Let x, y, (f) be real continuous functions on interval [a, b ], <p(t ) > 
0, t £ [a, b], and 



(42) 


Then holds 



(43) 


If, in addition, y(-) is a nondecreasing function defined on [a, b], we have 


x(t) <y(t)e^ mds . 


(44) 


Lemma 4 (@) Let x : [a,b] —> [0, +oo) be a real function and y(-) is a nonneg¬ 
ative, locally integrable function defined on [ a,b], a e (0,1), and there exists a 
constant L > 0 such that 



(45) 


Then there exists a constant C = C(a) such that 



(46) 


for all t G [a, b]. 
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Theorem 3 Under the assumptions given in Theorem [7] let u(t) and v(t) be the 
solutions of the Cauchy type problem (l20l) with different initial conditions. Then 


n ~ -go\e At 


\u(t) - v(t)\ < 


{it - a) a 1 + C(a)(t - a) la ! ), 0 < a < 1, 


I Ja) 


(t-aT 


(47) 


(p(t) [e a (t - a) a + 11 e At , n - 1 < a < n, n > 2 , 


n —1 


where ip(t) = C(a ) ^ \g k - g k \(t - af k g k = [ a D \“ k \e Ax u(t))\[ =a , and 


k=0 


it = [ U 7r* _ v-M0)]L, 

Similarly, under the assumptions given in Theorem |2] for the problem (l2ll) . there 
exists 


f|c 0 -'c k \e~ At {\ + C— —— 0 < a < 1, 

l«(0 - *01 < / 1 9 (48) 

I <A(0 (e “ (t - a) a + 1 Je A , n - 1 < a < n,n > 2, 


n— 1 


where f(t) = C(a ) ^ | c k -c k \(t-a) k , 


Ck 


k=0 


d k „ 

L5? (e “ ( '» 


, and c k = 


4(^v(0) 


Proof. Suppose that u(t) and v(t) are any two solutions of the Cauchy type problem 
(l20l) with different initial conditions. From the equivalent integral formulation (1221) . 
we have 


( n— 1 


e At \ u(t) - v(t) | < 


Zk* ■ **l 


(.t - a) c 


-k-l \ 


\k =0 

+ 


r(a - k) 


< 


f (t-sy- l e As \f(s,u(s))-f(s,v(s))\ds 

r(«) Ja 

( n ~\ f+ „\a—k— 1 N 

Z k* ■ 


(49) 


(t - af 


+ 


v *=o 

C/,p 

rc«j 


r(or - k) 


fit-sy- 1 * 

a 


e As \u(s) - v(s)\ds. 


For 0 < a < 1, using the generalized Gronwall’s inequality with weak singular 
kernel given in Lemma |H we get 


e At \u(t) - v(t )| <(Jg 0 ■ 
+ C 


(t - af~ l 
k(o') 


go ~ go 


1(5 - a) 


a -1 


T(a) 


(;t - sf- l ds. 


(50) 


which implies 
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\u(t) - v(f)| < |g 0 - S°|^((r -af 1 + C J^ (s- af \t - sf 'dsj 


T(a) 

_ I __ | e~ At 
- go - go 77—■ 

1 T(of) 

For n — 1 < a < n, n > 2, taking 
x(t) = e Al \u(t) - v(f)|, y(t) = 
in inequality (1421) . we have 


(t - af- 1 + C 


Jja ) 2 

r(2or) 


0 t-a ) 2ff - 1 ). 


(51) 


/«-1 


£ 1 * 


'k ~ gk\ 


l(t - a)‘ 


t-k -1 \ 


V A=0 


r(or - k) 


, (pit) = (t- s) 


a -1 


e At \uit) - v(f)| 
n_1 dt-af- k - ] 


2a—k -1 r (<*) 

e a 


^ y \gk - gk\ —7— + V IgA- - m - O) 

to na-k) ^ F (2 a-k) 

n— 1 (a 

< da) | gk -gk\it ~ af~ k ~ l [eM^it - af + l)e~ At , 


k=0 


(52) 


where C(or) = max I r , ‘ + r *} a] . 

Jt=0,l,-,n-l lr(Qr -*) r(2 a ~ k) > 

With the similar method, we can prove the stability results for the problem (I2T1) . □ 


3.3 generalized cauchy problems 


In this subsection, we consider the ODE with m- term Riemann-Liouville tempered 
fractional derivatives 


LD a t Mt) = fit,uit), a Df M uit), a Df^uit), ■ ■ ■ , a D a r uAm ~ l uit)\ n — \ < a < n, 
\[ a D“- k -\e Ax uit))\_ a = g k , k = 0,1,2, ■ • ■ ,n - 1; 

(53) 

and the ODE with m -term Caputo tempered fractional derivatives 


f aD^uit) = fit, uit), c X lM uit), c a D a ; M uit), 
d k 

= c k , k = 0,1,2, ■ ,n - 1, 

t=a 


C 1 1 / j\\ 1 

,D t n - 1 < a < n. 


’ a 




(54) 
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where A, > 0, the real value a £ (n - l,ri),n £ N + ,/ = 1,2, ••• ,m - 1, and 
0 < ai < cr 2 < • • • < cr m _i < a;, m > 2. Similar to the Cauchy problems (l20l) and 
(|2TI) . we have the following lemma. 


Lemma 5 Assume that B is an open set in R m and f : (a, b) x B —> R is a function 
such that f(t, u, u\,u 2 ,..., w m _0 e L([a, b])for cdl (n, «i, W 2 , • w m -i) s B. 

(1) If u(t ) £ L([a, b\) is the solution of the initial value problem (1531) if and only if 
u(t ) is the solution of the Volterra integral equation of the second kind 


u(f) 

„-Au+ „\ct-k -1 


X 


AV 


f (t - a)° 


£=0 

+ 


r(a - k ) 


J_ f e~ Mt ~ s) (t - s) a ~\f{s, u(s), a D‘f M i 

r(«) 


u{s\ a D a s ^ 2 u{s\..., a D\ 




u{s))ds. 

(55) 


(2) //"/(/, u, u i, « 2 , ..., u m - 1 ) /.S’ continous, then u(t) is the solution of the initial value 
problem (l54l) if and only ifu(t) is the solution of the Volterra integral equation of 
the second kind 


u{t) 


n —1 


Z 


c k 


e At (t - a) k 
Y(k + 1) 


+ 


Hor) 


f 


-A(t-s) 


(t-s)° 


~ l f(s, u(s),^D a ' Ji u(s), 


c r V y 2 ,A 2 / s 

a D s u(s),..„ 


C j~^m— 1 1 


u{s))ds. 


Suppose that f(t, u \, u 2 ,..., u m ) is a continuous function satisfying the Lipschitz type 
condition 


1/(6 u\, u 2 , -,u m ) - f(t,v i,v 2 , ...,v m )| < C Lipg ^ I Uj - Vj\ (56) 

7=1 

for all t £ [a, b ] and {u\,u 2 , ..., w m ), Oh, v 2 ,..., v m ) € fi, where the Lipschitz constant 
C Lipg does not depend on t £ [a, b\. 

Theorem 4 Let B be an open set in R m and f : [a, £>] x B —> R a function such 

that f(t, u) £ L([a, bf) for any u £ B and be Lipschitz continuous. 

(1) There exists a unique solution u{t) to the Cauchy type problem (1531) in the space 
L aA ([a,b]). 

(2) There exists a unique solution u(t) to the Cauchy type problem (l54l) in the space 
AC n [a,b]. 
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Proof. Similar to Theorems [H and [2l We begin our proof from the integral equations 
given in Lemma \5\ We only prove the Cauchy type problem (1531) . Let t\ belong to 
(a, b ) such that the inequality 


Upg Lr (a + *■ 


holds.The operator corresponding to (1551) takes the form 


(Tu)(t) = u Q (t) + —— I e 

r(a) 


— f, 

(a) 1 




u(s), a D a s u l u(s ),..., a D a r'^-'u(s))d.s, 

(57) 


where 


n— 1 

«o(0 = git 


i=0 


- a) a ~ k ~ ] 
T(or - k) 


(58) 


From the Lipschitz condition (1561) it directly follows that 


u(t), ( a D^um, ..., ( a D a r'^-'um 
- v(t), ( a DX M vm, ..., { a D a r uAm -'v)m 

m— 1 


y*l»^l, 


V^Vn-T ,A m —i 


< C I 

— Lipg a L t 
m— 1 

Z(T 


Y^aD^iu-v) 

7=0 


(0 


< CV.,: 




a—or/,/il r Q;i,A r^ai-A;, , \ , xx i i 

J I//’ a D t \u(t)-v{t)) . 


7=0 


Furthermore, using the composition formula (fl2l) . we have 


- jXfit, v(f), a D? M v(t ),..., 

m—1 

— Cupg 




V(f))l 


m-1 

Z(»L“' 


7=0 

— 1 e~ At (t - a) aj ~ kj ~ l 


Z 

k:=0 


n a J - k >) 


(w - v) 

\ a D a t rkr \e A \u{t) - v(0))| f= J )(0, 


where is the smallest integer larger than or equal to aj. 
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From the given initial conditions, it can be checked that a D® ] kj ' (e A '(u(t)-v(t)))\ = = 
0. Then 


Ur J fit, U(t), aD? lM u(t), ..., a D a r" Xm -'u(t)) 

- a r; A f(t, v{t), a D? M v(t ),..., a D a r'^-'vm 

m-i 

<c Lipg 2( a /r a "i(«(o-v(o)|), 

7=0 


(59) 


for any t e [a, b]. Taking t = t\ in above the formula and applying (IA.2I) . we get 


\\TU\ - T M 2 ||l([o,/ 1 ]) ^ Cu P gK\\U\(t) - U2{t)\\L([a,t l }), K - Cupg 


in— 1 


z 


(fi - ay~ a > 

r(a - ctj + 1) 


It follows that there exists a unique solution u*{t ) to Eq. (l55l) in L([a, ti]). This 
solution is obtained as a limit of the convergent sequence ( T j u* 0 )(t ) = «,(/), and 
holds 


lim || T J u* 0 - u*\\u[ aM) = 0, (60) 

7 -*°° 


i.e., 


lim || uj - u*\\ L([aM]) = 0. (61) 

j->°° 


With the same fashion of proving Theorem [H we can show that there exists an 
unique solution u(t) e L a ’ A ([a,b ]) to Eq. (1551) . In addition, 


I a D°' A Uj - a D a t ’ A u\\ L qa^) 

= \\f(t,Uj-U a D“' M Uj-i, .. 
- fit, u, a D a t lM u, ..., a D" 


a D a r uK -' UjH ) 

W)|| L([a,b]) 


< K\\Uj-l - U\\u[a,b\) 0. j -» 00 , 


which implies that a Df A u(t) e L{[ci, b ]). 


□ 


4 Numerical algorithm for the tempered fractional ODE 


The well-posedness of the tempered fractional ODE has been carefully discussed 
in the above sections. Usually, it is hard to find the analytical solutions of the tem¬ 
pered fractional ODE, especially for the nonlinear case. Efficient numerical algo¬ 
rithm naturally becomes an urgent topic for this type of equation. Now, we extend 
the so-called Jacobi-predictor-corrector algorithm |44| to the the tempered frac¬ 
tional ODE; and its striking benefits are still kept, including having any desired 
convergence orders and the linearly increasing computational cost with the time t. 
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4.1 Jacobi-predictor-corrector algorithm 


The equations (l22l) and (l28l) can be rewritten as the following Volterra integral 
equation 


p, e~ M r 

u(t) = e At y a k (t) + —— 

U r ( ff ) 


(t - s) a l g(s, u(s))ds. 


(62) 


where a k {t) = c k 777777 for tempered Caputo derivative, and a k (t) = gk ' p (f _ k) for 
tempered Riemann-Liouville derivative, and g(t, u(t )) = e At f(t, u(t)). 


Firstly, we transform the original equation into 


u( , ) = e-Xa M+ pf 


(t - s) a l g(s, u(s))ds 


k =0 
n —1 


= e 


-At 


a k(t) 


+ 


k=0 


e (t-a > 


) J' (1 -zf l g(z,u(z))dz, 


(63) 


where 


g(z, u(z )) = g 


a t + a jt 

—z + —— ,u\ 


a t + a 

—z + —— 


- 1 <z< 1; 


N (t-a t + a\ 

u(z) = u (—2~ z + ~Y~ J’ - 1 ^z<l. 

Using (N + l)-point Jacobi-Gauss-Lobatto quadrature to approximate the integral 
in (1631) yields 


rl— I —At . rv IV 

^ ^ at / ^ ^ \ u ^ 

uui - ) at U) + —-(—— > «Jj g(Zj,u(Zj)), (64) 

u r(a)l 2 > U 

where we choose co(z) = (1 -z)" _1 (l + z)° as the weight function; {Zj}J =0 and {a>j}J =0 
are the (At + l)-degree Jacobi-Gauss-Lobatto nodes and their corresponding Jacobi 
weights in the reference interval [-1,1], respectively; see, e.g., If38ll . 

Now we turn to describe the computational scheme for Eq. (l63l) . For this purpose, 
we define a grid in the interval [ a , b ] with M + 1 equidistant nodes tj, given by 

tj = jt + a, j = 0, • • • , M, (65) 

where r = (b-a)/M is the stepsize. Suppose that we have got the numerical values 
of u{t ) at t 0 , t\, • • • , t n , which are denoted as w 0 , Mi, • • • , u n ; now we are going to 
compute the value of u{t) at t n+ 1 , i.e., u n+ \. From Eq. (l64l) . we have 
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w(tn+l) 


n— 1 


~e ' K+ 1 ^ j a k (t n+ 1) + 

k =0 
/2— 1 

= e~ /lf,,+1 ^ a*(4+i) + 


1 tfl+\ 


N 


k=0 


r (a) 

\ 2 

7 7=1 

1 / 

'4+i - a' 

a N 

) E 

7=0 

r (a)\ 

, 2 , 


^ ^ Sn+X^Zj-, U n+ l(Zj)) 


u)je" + '- A(Zj l) f n+1 (zj,u n+ i(Zj)),(66 ) 


where 


r / ~ /./4+1 ® 4+1 ^ 

/n+l(z, U n+ i{z)) = /(--- Z + - r -, U 


tf2+ 1 ^ 1 I 

- ^ - £ + 


■F))' 


1 < Z < 1, 


*4+1 (z) ** 


4+1 — O 4+1 + ^ 

-r- Z + --- 


, - 1 < Z < 1. 


To compute the second summation term of (l 66 l) . we need to evaluate the value of 
/ at the point t n+] due to f„+\(zN, * 4 +ifev)) = /( 4 +i, w( 4 +i)). It can be numerically 
approximated by using the piecewise linear interpolation of the term / H14I13I10L 
Here, we do it using the technique given in our previous work P4l . More con¬ 
cretely, we interpolate the function / by the known “neighborhood” points of t n+ 1 . 
For the other values of f n+ \(zj, w«+i(z/)), 0 < j < N - 1, we can also obtain them 
based on the interpolation of / at the time nodes located in the “neighborhood” 
points of Zj (should be (1 + Zj)t n + 1/2 as to variable t). Denote Nj as the number of 
time nodes used for the interpolation. In practical application, we use the improved 
predictor-corrector formulas given in IfTOTl to get the values at to, t\, ■ ■ ■ , t Nl -\ as the 
known ‘initial’ values. For the criterion of choosing the “neighborhood” points, 
refer to ||44| . 

Collecting the above analysis, we get the predictor-corrector formulas of (l62l) as 


Wn+1 


n—l 


e Atn+< Y a k (t n+] ) 


k=0 


+- 


1 (t n +i ** 


Ha) 


f N -1 

(E 

j =o 


coje 


~ A{Zi 1} f n +l(Zj,U n+1 (Zj)) + CU/v/(4+1,^ +1 )),(67) 


and 


n -1 


Kl = e~ Atn+X Yj a ^ tn+l) + fj~)( r " + 2 a ) Y w j 6 ^ 1] fn + l(.Zp &n+l(Zj)), 


k =0 


N 


7=0 


(68) 
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where {f£ +l j}*f =0 in (1681) means that all the values of f n +\ at the Jacobi-Gauss- 
Lobatto nodes are got by using the interpolations based on the values of {/(Y„ m,)}" = 0 ; 
whereas {f n +\j} ! Jj ( ' in (f67l) are obtained by using the interpolations based on the val¬ 
ues of {f(tj,Uj)}'j =0 and f(t n+ i, w^ +1 ). From the computational scheme (I67l)-(l68l). it 
can be clearly seen that the computational cost linearly increase with n (or time t). 

With the similar methods given in ||44j|. we can get the following estimates for the 
Volterra integral system (l62l ) . 

Theorem 5 If git, u(t)) is Lipschitz continuous with respect to the second variable, 
and has the form 

m 

git, u(t)) = ^ f k w k (t) + S, 

k= 1 

where w k (t), 1 < Ni < p\ < ... < p,„ are sufficiently smooth and m can be +oo, 6 is 
a constant, then there exists a constant C being independent ofn, r, N, such that 

max \u(t n+ \) - u n+ 1 | < Cr N ‘. 

l<n+l<M 1 1 


This theorem shows that the scheme (I67l)-(l68l) potentially have any desired conver¬ 
gence order by adjusting the number of interpolation points Nj. 


4.2 Numerical test 


In this subsection, we consider two simple numerical examples to show the numer¬ 
ical errors and convergence orders of the Jacobi-predictor-corrector method. The 
two examples are Caputo tempered ODEs; solving the Riemann-Liouville tem¬ 
pered ODEs can be done in the same way, so is omitted here. 

Example 1 Consider the Caputo tempered factional initial value problem 

o DfMt) = e~*( r( ^ ) g) f 8 ~ g + t* + 9 -C + ^T(a + 1)) - u(t). (69) 

The initial values are chosen as u{t)| r=0 = 0 and [j t (e M u(t ))\| = 0 for 1 < a < 2, 

and as u(t )\ t= o = 0 for 0 < a < 1. Using the formula 

< 70) 

it can be checked that the exact solution of this initial value problem is 

u(t ) = e~ A, (t* + 
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Table 1 

Maximum errors and convergence orders of Example Q] solved by the scheme (l67l)- (l68l) 
with T = l,N = 20, Ni = 7, and a = 0.5. 



T = 0 

A = 2 

A = 6 

T 

error 

order 

error 

order 

error 

order 

1/10 

1.5207e-004 


2.3516e-005 


1.4300e-006 


1/20 

4.6202e-007 

8.3626 

1.4040e-007 

7.3879 

3.3507e-008 

5.4154 

1/40 

1.6877e-009 

8.0967 

6.3106e-010 

7.7976 

2.7846e-010 

6.9109 

1/80 

8.1135e-012 

7.7005 

2.5491e-012 

7.9517 

1.4371e-012 

7.5982 

1/160 

3.5305e-014 

7.8443 

1.2794e-014 

7.6383 

7.0913e-015 

7.6629 


Table 2 

Maximum errors and convergence orders of Example Q] solved by the scheme (I67l)-(l68l) 
with T = 1,1V = 20,Ni = 6, and a = 1.0. 



4 = 0 

4 = 2 

4 = 6 

T 

error 

order 

error 

order 

error 

order 

1/10 

8.1108e-005 


1.2528e-005 


1.1365e-006 


1/20 

7.8788e-007 

6.6857 

1.5673e-007 

6.3207 

2.3299e-008 

5.6082 

1/40 

1.2817e-008 

5.9418 

2.1909e-009 

6.1606 

3.2657e-010 

6.1567 

1/80 

2.2418e-010 

5.8373 

3.4124e-011 

6.0046 

4.4768e-012 

6.1888 

1/160 

3.6193e-012 

5.9528 

5.3461e-013 

5.9962 

6.7955e-014 

6.0417 


Table 3 

Maximum errors and convergence orders of Example Q] solved by the scheme (|67T)-(|68T) 
with T = 1,1V = 20,Nj = 6, and a = 1.5. 



4 = 0 

4 = 2 

4 = 6 

T 

error 

order 

error 

order 

error 

order 

1/10 

6.6386e-005 


9.6009e-006 


8.5068e-007 


1/20 

9.2847e-007 

6.1599 

1.4297e-007 

6.0694 

1.9943e-008 

5.4147 

1/40 

1.5767e-008 

5.8799 

2.1338e-009 

6.0661 

3.0437e-010 

6.0339 

1/80 

2.3505e-010 

6.0678 

3.5138e-011 

5.9242 

3.8203e-012 

6.3159 

1/160 

3.8498e-012 

5.9320 

5.3434e-013 

6.0391 

6.7433e-014 

5.8241 


In our numerical algorithm, the values of u(t) at points {/ 7 }^ 0 are calculated by the 
improved Adam’s methods IITOll . The numerical results are reported in Tables dB 
From Tables QB we can see that the convergence orders are in good agreement 
with the theory presented in Theorem [5} 
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Example 2 In this example, we examine the following initial value problem 


c 0 Df A u(t) = -pu(t), p>0. (71) 

The initial values are given as e A, u{t)| f=0 = 1 and \j t (e At u(f ))\| = 0 for a e (1,2), 

and as e A, u{t) | f= o = 1 for a e (0,1). Using the Laplace transform presented in 
Appendix B, we have 

(.s + A)°7l(s) -(s + Af- 1 = -p u(s). (72) 


Then 


u{s) 


Employing the Laplace transform 
function l \33\] 


(s + A) 


a -1 


(S + A)° 


(73) 


involving the derivative of the Mittag-Leffter 


Z{t p -\E afi (~af))} 


s a ~P 
s a + a' 


Re(s ) > \a\ l,a , 


(74) 


we can check that the exact solution of this initial value problem is 


u{t) = e At E a f-p f). 


Here the generalized Mittag-Leffter function E a ^(-) is given by l \J3\l 


E a p{z) = Yj 

k =o 


Y(ak+/3y 


Re(a) > 0. 


(75) 


In this example, the solution u{t) does not have a bounded first (second) derivative 
at the initial time / = 0 as 0 < a < 1(1 < a < 2). To improve the convergence order, 
we employ the technique given in our previous work 1.441 . We separately solve the 
equation in subintervals [0, T 0 \ and [ 7 0 , T\ of the interval [0, 7’]. More specifically, 
we modify the formula (1631) as 


e~ At r T ° 

u(t) = e~ Al V a k {t) + —— (t - sf^g{s, 

fli E(a) Jo 


u(s))ds 


k=0 
-At nt 


+ - 


/< 

Jt 0 


= e 


r(or) 

n— 1 

At Y a k(0 

k =0 
-At 


{t - s) a g(s, u(s))ds 


-At r T « 


+ 


r(a) 


f 


(t - s) a g(s, u(s))ds 


p—At / f _ p \ a /"•! 

+ rw(^) 
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Here, we suppose that the smoothness of g is weaker on the subinterval [0, T () \ and 
sufficiently smooth on \T (h T\. For the integral on the subinterval [0, 7 q], the Gauss- 
Lobatto quadrature with the weight function w(s) = 1 is used; and for the one on 
[T 0 , T ], we compute it as ([64b . i.e., 


n- i i N 

X(t) * ^ a k (t) + —- ^ d)j(t - Sj) a ~ l g(sj, u(Sj )) 


k=0 


r(a)\ 2 


7=0 

N 




7=0 


(76) 


where /V, and {v ; }^ =0 correspond to the number of, the weights of, and the 

values of the Gauss-Lobatto nodes with the weight oj(s) = 1 in the interval [0, 7o], 
respectively. The values of {g(sj, v(Sj))}^ =0 can be computed as in the starting pro¬ 
cedure. Since g and x are continuous in the interval [0, T 0 ], by the theory of Gauss 
quadrature llBTll and the analysis above, we can see that if N is a big number then 
the accuracy of the total error can still be remained. The numerical results are re¬ 
ported in Tables 0 and [5} And it can be seen that the desired numerical accuracy is 
obtained. 

Table 4 

Maximum errors and convergence orders of Example 0 solved by the scheme (l76l) with 
T = 1.1,1V = 20,1V - 40, Ni = 2, T 0 = 0.1, n = 1, and A = 5. 



a - 0.2 

a - 0.9 

a = 1.8 

T 

error 

order 

error 

order 

error 

order 

1/20 

5.4805e-004 


1.9043e-005 


2.146le-006 


1/40 

1.8749e-004 

1.5475 

4.3478e-006 

2.1309 

5.6685e-007 

1.9207 

1/80 

5.0838e-005 

1.8828 

1.085 le-006 

2.0025 

1.5416e-007 

1.8786 

1/160 

1.3492e-005 

1.9138 

3.1549e-007 

1.7821 

4.0386e-008 

1.9324 


5 Concluding remarks 


Currently, it is widely recognized that fractional calculus is a powerful tool in de¬ 
scribing anomalous diffusion. Because of the bounded physical space and the finite 
lifetime of living particles, in the CTRW model, sometimes it is necessary to trun¬ 
cate (temper) the measures of jump length and waiting time with divergent second 
and first moments, respectively. Exponential tempering offers technique advantages 
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Table 5 

Maximum errors and convergence orders of Example |2] solved by the scheme (l76l) with 
T = 1.1,1V = 20,1V = 40,1V/ - 2,Tq - 0.1,yu = 1, and A - 10. 



a = 0.2 

a = 0.9 

a = 1.8 

T 

error 

order 

error 

order 

error 

order 

1/20 

2.0162e-004 


7.0054e-006 


4.0825e-007 


1/40 

8.8563e-005 

1.1868 

1.6897e-006 

2.0517 

1.0286e-007 

1.9887 

1/80 

2.721 le-005 

1.7025 

4.1757e-007 

2.0167 

2.7730e-008 

1.8912 

1/160 

7.4508e-006 

1.8688 

1.1169e-007 

1.9025 

7.3208e-009 

1.9214 


of remaining the infinitely divisible Levy process and that the transition densities 
can be computed at any scale. Then the Fokker-Planck equation of the correspond¬ 
ing CTRW model has the time tempered fractional derivative (tempered waiting 
time) and/or the space tempered fractional derivative (tempered jump length). This 
paper focus on discussing the properties of the time tempered fractional deriva¬ 
tives as well as the well-posedness and numerical algorithm for the time tempered 
evolution equation, i.e., the tempered fractional ODEs. The proposed so-called 
Jacobi-predictor-corrector algorithm shows its powerfulness/advantages in solving 
the tempered fractional ODEs, including the one of easily getting any desired con¬ 
vergence orders by simply changing the parameter of the number of the interpolat¬ 
ing points and the other one of linearly increasing computational cost with time t 
rather than quadratically increasing more often happened for numerically solving 
fractional evolution equation. 
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Appendix A: Properties of tempered fractional integral and derivatives 


For the Riemann-Liouville tempered fractional integral, we have the following 
properties. 
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Property 1 Ifu(t) e AC n [a, b ], cr > 0, and d e R, then for all t > a, 


J?Mt) = X 

£=i 


e~ A, (t - aY +k ~ l 

d k ] (e At u(t)) 

- 

r(cr + k) 

1 

t=a- 




+ 


r(cr + n ) 


f 


d" 

(t - 5) cr+ ' ,_1 ^ 


— (e As u(s))ds. 
ds n 


Proof. If u(t) has continuous derivative for t > a. then using integration by parts to 
CD, there exists 


i- [\,-,r-‘e-u(sVsA^m + 

r(cr) J a 


r + 1 ) Ja ds 


e As u(s)ds; 


r(cr + 1) T(cr 

and if the function has n continuous derivatives,then integrating by parts, we get 
J 0 - s) a ~ 1 e As u(s)ds = ^ 


k= 

+ 


1 

+ 

1 

\d k -\e At u{t)) 

- 

r(cr + k ) 

dt k ~ l 

t-a- 


no- 


-- f (t - sf +n 1 ^—(e As u(s))ds. 

+ n) J o ds n 


(A.l) 


Multiplying both sides of (IA.1I) by function e At , the desired result is obtained. □ 


In addition, integrating by parts leads to 

u{s)ds + I" e~ cr( - x ~ s ' > du(s), 

J a 

which implies that if u(t ) has continuous derivatives on finite domain [a,b], and 
cr > 0, then for all t > a, 

lim XMt) = u(t). 

cr-»0 




cr,A 


u{t) = e ^u^a) + cr 


/' 


,-<T(t-S ) 


Lemma 6 For a function u{t) e L([a, b ]), cr > 0, d e R, we have 

||_^_ J e ^ s)(t _ sr -i u(s)ds ^^ < M\\u\\ U[aM) , 


where M 


■ 

r(cr+l) - 


(A.2) 
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Proof. By simple calculation, we have 



e ~A(ts) 


(,-sr-'u(s)ds\f 


?- Ait ~ s \t - sy- 1 


a a 
b n b 


1 r r 
< 1^) 1 Js 




< 


T(cr) 

(b - aY 
r (<t + 1) 


INI L([a,b\)- 


\u(s)\dsdt 

dt\u(s)\ds 


□ 


Proposition3 (Semigroup properties) Let u(t) £ L(\a, b\) and cri,<x 2 > 0, A £ R. 
Then for all t > a, 


rO-i.A 

a 1 t 


\j?Mt)\ 


i: T,+<T2j U (t) 


at 


iyXiyMm. 


at 


Proof. Recalling the semigroup property of Riemann-Liouville fractional integral 
11331 If Lpuit) = Ip +(T2 u(t), we get the following semigroup property of tempered 
fractional integral 

jyXiyMm = e- M a iy\e%iyMm 

= e- At jy[e A \e- At a iy(e A, um\ 

= e- At a iy[ a iy(e A, u(t))\ 

= e - At a ry +CT2 (e At u(t)) 

= a iy +CT2 Mt). 


□ 

Similar to the fractional calculus, the tempered fractional calculus is a linear oper¬ 
ation. 

Proposition 4 (Linearity properties) Let u(t) £ L([a, b\ ). Then for all 6,p,A £ R: 

( 1) for <j > 0, holds 


jy\5 u(t) + p U (t )] = d jy\u(t)} + p jy A [um 


(2) for a £ (n - 1 ,n), holds 

a Df A [6 u(t) + p u(t )] = <5 a D a /[u(t)] + p a Df\u(t)l 
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(3) for a £ (n - 1 ,n), holds 


c 

a 


Df\6 u(t ) + p u(t)] = 6 c a D a t \u(t )] + p c a Df\u(t)l 


Proof. The linearity of fractional integral and derivatives follows directly from the 
corresponding definitions. We omit the details here. □ 


Propositions Let u(t ) £ C"[a,b] and a £ (n - \ ,n). Then for all t £ {a,b\ there 
exists 


c 

a 


D a f\um 


t—Cl 


= 0 . 


Proof After simple argument, we have 

”- At r ‘ 1 d n e As u(s ) 


< 


e 

Rh 


I c .d"; "«| 

where M = max r€ [ a ^ 


- At r r l 

1 (t - sy-" +1 


ds n 


ds 


< 


Me~ At {t - a) ,, -‘ 
T{n - a + 1) 


(A.3) 


I d n e x, u(t) 
I dt n 


. From the above analysis, we get the desired result. □ 


Appendix B: Laplace transforms of the tempered fractional calculus 


In this subsection, we discuss the Laplace transforms of the tempered fractional 
calculus. Define the Laplace transform of a function u(t) and its inverse as 071 


£{u(t); .v} = u(s) 


=JT 


e st u{t)dt. 


(B.4) 


and 


1 /^C0+lOO 

£~ l {u(s); t} = u(t) =—. I e sr u(s)ds, Co = Refs) > 0, i 2 = - 1. (B.5) 

co—ioo 

We start with the Laplace transform of the Riemann-Liouville tempered fractional 
integral of order cr. 


Proposition 6 The Laplace transform of the Riemann-Liouville tempered fractional 
integral is given by 


£( 0 If J u(t)) = (A + s) <r u(s). 


(B.6) 


Proof. First, we rewrite the Riemann-Liouville tempered fractional integral as the 
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form of convolution 


oCVo 


=-Lr, 

Her) J„ 


' s) (f - 5)°" l u(s)ds = 


r(m) 


M(f). 


In view of the Laplace transform of the convolution 071 

X{w(0 * v(f); s} = u(sjv(s). 


we have 


x { 0 itMty, s} = £ 


e -At t <r-l \ 

s\ £{u(t); s}. 
Her) J 


Recalling the Laplace transform £{e At f 1 ;5} = r(cr)(/l+5) <r , we have the Laplace 
transform of Riemann-Liouville tempered fractional integral 


£{ 0 If A u(t); 5 } = (A + s) cr u(s). 


(B.7) 


□ 

Next, we turn to consider the Laplace transform of tempered fractional derivative. 

Proposition 7 The Laplace transform of the Riemann-Liouville tempered fractional 
derivative is given by 


£{ 0 D a /u(t)-s} 


is + A) cr u(s) - + A) k 


k =0 


o D a - k -\e xt 


u(t))\ t= 


(B.8) 


The Laplace transform of the Caputo tempered fractional derivative is given by 

n- 1 _ jk 

t =0 


k=0 


£{%D a f A u(t)\ s} = (s + A) ir u(s ) - ^0 + A) a k 1 ^(.e At u(t))[ 


(B.9) 


Proof. Using the properties of Riemann-Liouville tempered fractional calculus, we 
may rewrite the Riemann-Liouville tempered fractional derivative as 

oDf A u(t) = e~ A, ^p(v(t)), n - 1 < a < n, 

where v(t) denotes 


v(t) = 


Tin 


— f 

- a) Jo 


^^o/rvvo). 


«- .v) 


Furthermore, using the formula of Laplace transform of an integer order derivative 


n -1 




k =0 


d^ 

i? v<,) 


n— 1 


t =0 


v—i 1.1 d" ^ 1 

= ."7W - ^ —v(l) 


k =0 


/=0 
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Combing the first translation Theorem 071 


£{e At u(t)', s} = u(A + s), A £ R, Re{s ) > A, 


(B.10) 


we have 


£ \ s \ = (s+Ams + A) ~ Tj (s+A) k 


n— 1 


in—k —1 


dt 

which implies 


*=o 


n— 1 


d t n 


r An-k-l 


t= 0 


£{oD?Mty, S} = (s + T)"T(s + A)~ 


k =o 


i=0 


By applying the convolution theorem once again, we get 


£{v(ty,s} =£{ 0 I'r(e At um 


-£ 


+n-a—l 


; s\ £{(e At u(t))\ 5 } 


Y{n - a) 

=s- {n - ay n(s - a). 

Combining above formulae, we have 

V(s + A) = (s + A)~ (n ~ ar u(s). 


(B. 11) 


(B.12) 


(B. 13) 


In addition, from the definition of the Riemann-Liouville fractional integral it fol¬ 
lows that 


in-k-l 


dt 


<n—k— 1 


d ^' <0 = jF^ d,rie " U(,)) = 

Substituting (IB.13I) and (IB. 141) into (IB. 121) . we have 

n —1 

£{ 0 D c ;Mt); S} = (s + A) or u(s) - + A) k ()Dy (k ~ a+]> (e Ax u(t))\ i=Q . (B.15) 


(B.14) 


k=0 


To establish the Laplace transform for the Caputo tempered fractional derivative, 
let us write it in the form 

c 0 Dfum =e WAd)) = 0(0). 

V(,) = l> (,)) ' 

w(t) =e At u(t). 

Applying the Riemann-Liouville fractional integral and the first translation Theo¬ 
rem (IB. 101) . we then have that 

£{ c Q D a t ' A u{ty s} = (s + Ay {n - a) Y( S + A), (B.16) 
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where 


n —1 

v(s + A) =(s + Afw(s + A) - J^(s + Af~ k ~ l 

k=0 


d k 

d? w{,) 


f=0 


(B.17) 


Noticing that £{w(t); s} = u(s - A), and inserting (IB.171) into (IB.161) . we arrive at 
the Laplace transform of the Caputo tempered fractional derivative (IB.9I) . □ 


It is easy to see that if taking A = 0 in (IB.81 ). we have the Laplace transform of 
Riemann-Liouville fractional derivative If33l 


n— 1 

£{ 0 D°u(t); 5 } = s cr u(s) - ^ s>c 

k=0 

And if taking A = 0 in (IB.9D . we can get the Laplace transform of Caputo fractional 
derivative 031 


J) 


-(k-a+1 ) 




(B. 18) 


n -1 


X{qZ)“m(0; 5 } = s ,r u(s) - ^ 


i=0 


k-1 

d k u(t) 

- 


dt k 

t=0 


n— 1 < a < n. 


(B.19) 


From the Laplace transform of tempered fractional derivatives, we observe that dif¬ 
ferent initial value conditions are needed for fractional differential equations with 
different fractional derivatives. From (IB.9I) . it can be noted that the Laplace trans¬ 
form for the Caputo tempered fractional derivative involves the values of the func¬ 
tion u{t) and its derivatives at the lower terminal t = 0, which are easily specified in 
physical. So the Caputo type fractional derivatives are more popularly used in time 
direction Il33l . The Laplace transform for the variants of the Riemann-Liouville 
tempered fractional derivatives are given as 


(.A + s) a u(s) - [q/, 1 a (e At u(t))\ r=0 ] - A“u(s), 0 < a < 1, 


M'VO) 


(d + sflAis) - aA a l (su(s) - u{ 0)) - A LY u{s) - [ 0 D a t 1 (e At u(t))\ t=0 \ 
-(A + s)[ 0 I?- a (e A *u(0)l t _ o \, l<a<2. 

(B.20) 
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